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ABSTRACT 

We briefly describe a library that may be used with any e + e~ — > T + T~(n^) Monte Carlo 
program to account for the effects of anomalous tt'j couplings. The implementation of 
this library in KORALZ version 4.04 is discussed. 
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1 Introduction 



Radiative r pair production is of great interest, as it is sensitive to anomalous electro- 
magnetic couplings of the r. With the sensitivity afforded by the LEP experiments, this 
provides an opportunity to search for new physics phenomena [1-4]. Any meaningful 
interpretation of the experimental data requires a Monte Carlo simulation in which Stan- 
dard Model predictions may be augmented by the contributions from possible anomalous 
couplings. 

Since the LEP collaborations are entering their final years of operation it is a good 
time to document the programs that were actually used in data analyses. In this paper 
we describe a library that has been used to calculate anomalous contributions to TT7 
couplings [[|. The library is based on the work described in Q and can be used with any 
e + e~ — > T + T~(n , -f) Monte Carlo program and, after minor adaptation, withpp — > Z/ r y+...; 
Z/ry — > T + T^{n r )) or ep — > Z/7 + Z/7 — > T + T~(n r y) programs as well. 

In the present paper, we will discuss the interface of our library to KORALZ version 
4.04, which is described in detail in || |J. The fortran code of the library is archived 
together with KORALZ ||, in the same tree of directories. Let us note that in the 
future, KORALZ will be replaced by a new program, KK2f 0, which is based on a more 
powerful exponentiation at the spin amplitude level; implementation of our library will 
be straightforward for that program as well. 



2 Calculation of anomalous couplings 

To evaluate the effects of anomalous electromagnetic couplings on radiative r pair produc- 
tion, a tree-level calculation of the squared matrix element for the process e + e~ — > t + t~7 
has been carried out @], including contributions from the anomalous magnetic dipole mo- 
ment at q 2 = 0, -^2(0), and the electric dipole moment F 3 (0). This calculation is included 
in our library. When activated, it uses the 4-momenta of the leptons and the photon 
generated by the host program to compute a weight, w, for each event according to 



\M 



I ano / .« \ 

W = TT1 il- U) 



A4 ano is the matrix element for F 2 (0) 7^ and/or F 3 (0) 7^ 0, and A^sm is the matrix 
element for F 2 (0) = F 3 (0) = 0. 

As this calculation is performed at 0(a), the case of multiple bremsstrahlung requires 
special treatment. In this case, a reduction procedure is first applied in which all photons, 
except the one with the greatest momentum transverse to a lepton, or p?, are incorporated 
into the 4-momenta of effective initial- or final-state leptons. The 4-momenta of the 
photon with greatest pt and the effective leptons are then used to compute the weight. 
Cross-checks of the calculation against an independent and slightly simplified analytical 
calculation || as well as checks of the validity of the reduction procedure are described 
in [[J. The results of the calculation have been used in the measurement of anomalous 
electromagnetic moments of the r described in H. 
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3 Flags to control anomalous couplings in KORALZ 



In KORALZ version 4.04, the calculation in our library is activated by setting the card 
IFKALIN=2. This is transmitted from the main program via the KORALZ input parameter 
NPAR(15). Additional input parameters are set in the routine kzphynew(XPAR,NPAR), 
although there are currently no connections to the KORALZ matrix input parameters 
XPAR and NPAR. Table [I] summarizes the functions of these input parameters. 



Parameter 


Description 


Default 


IFL1 


Compute weights for ^(0) if IFL1 = 1 


1 


IFL2 


Compute weights for F 3 (0) if IFL2 = 1 





ISFL 


For ISFL = — 1, compute only terms with anomalous contribu- 
tions. For ISFL = 0, include all terms (anomalous, Standard 
Model, all interference). For ISFL = 1, use the approximation 
of ref. [p|] 





IRECSOFT 


Generate only events with photon(s) if IRECSOFT = 1 





EMINACT 


Minimum sum of all photon energies required to calculate 
anomalous weights 


17 GeV 


EMAXACT 


Maximum sum of all photon energies allowed to calculate 
anomalous weights 


1000 GeV 


PTACT 


Minimum sum of all photon momenta transverse to the beam 
direction required to calculate anomalous weights 


2 GeV 



Table 1: Input parameters to control the calculation of weights for anomalous electro- 
magnetic moments. 



In order to provide the user with enough information to retrieve w for a given event 
for any F 2 (0) or F 3 (0), we take advantage of the fact that, for each event, one may write 
w as a quadratic function of the anomalous couplings: 

w = aFf (0) + /3F 2 (0) + 7 F 3 2 (0) + 5F 3 (0) + e . (2) 

When the calculation of w is completed, the 5 weight parameters a, [3, 7, 8 and e are 
stored in the common block common /kalinout/ wtkal(6), with the assignments shown 
in Table 0. The user is then free to calculate w for whatever combination of F 2 (0) 
and -^3(0) is desired. In practice we set e = 1, since anomalous terms must vanish for 
F 2 (0) = -F 3 (0) = 0. We also set 5 = 0, as the interference between Standard Model and 
anomalous amplitudes vanishes in the case of radiation from an electric dipole moment. 
These short cuts save substantial CPU time. 
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Common block entry 


Weight parameter 


wtkal(l) 


not used here 


wtkal(2) 


e 


wtakl(3) 


a 


wtkal (4) 


P 


wtkal(5) 


7 


wtkal (6) 


5 



Table 2: Correspondence between entries in kalinout common block entries and weight 
parameters of eq. (H). 



The code for calculation of the weight w is placed in the directory korz_new/ttglib 
in the files ttg.f and ttgface.f. 

4 Demonstration programs 

The demonstration program DEM02 . f for the run of KORALZ when our library is activated 
can be found in the directory korz_new/f ebruary and the output DEM02.out in the 
directory korz_new/f ebruary/prodl. The Standard Model DEMO.f for KORALZ and 
its output DEMO . out are also included in the directories mentioned above. All these files, 
as well as the library itself, are archived together with KORALZ ||. The maximum 
centre-of-mass energy allowed for the runs with anomalous TT7 couplings is 200 GeV. 
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